Hopping transport in systems of finite thickness or length 
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Variable-range hopping transport along short one-dimensional wires and across the shortest di- 
mension of thin three-dimensional films and narrow two-dimensional ribbons is studied theoretically. 
Geometric and transport characteristics of the hopping resistor network are shown to depend on 
temperature T and the dimensionality of the system. In two and three dimensions the usual Mott law 
applies at high T where the correlation length of the network is smaller than the sample thickness. 
As T decreases, the network breaks into sparse filamentary paths, while the Mott law changes to a 
different T-dependence, which is derived using the percolation theory methods. In one dimension, 
deviations from the Mott law are known to exist at all temperatures because of rare fluctuations. 
The evolution of such fluctuations from highly-resistive "breaks" at high T to highly-conducting 
"shorts" at low T is elucidated. 

PACS numbers: 72.20.Ee, 73.63.Nm 



I. INTRODUCTION 

Variable-range hopping (VRH) is the mechanism of 
low-temperature transport common in systems with elec- 
tron states localized by disorder. Typically, the VRH 
conductivity (t{T) obeys the Mott law, 



lnCT(r) = - 



CdTn 



i/(d+i) 



Tn- 



1 



(1) 



where d is the space dimension, a is the localization 
length, g is the density of states, and Cd is a numeri- 
cal coefficient. Temperature dependence stronger than 
the Mott law can arise due to electron interactions that 
deplete g near the Fermi level.l^ On the other hand, T- 
dependences weaker than the Mott law have also been 
observed. For example, a power-law scaling. 



a{T) - T", 



(2) 



has been measurecpEl jn transport across disordered thin 
films. Such a behavior is thought to originate from hop- 
ping along special highly conductive chains of sites. SHZl 
The exponent a scales with the number of hops in the 
chain. We call this transport mechanism rare-chain hop- 
ping (RCH). The RCH has been often discussed in the 
context of magnetic tunneling junctions, see a recent ex- 
ample in Ref . HI 

In the preceding papei^l we have shown that the RCH 
can also determine the conduction measured in an ensem- 
ble of short one-dimensional (ID) wires connected in par- 
allel. In our theory the exponent a depends on the ratio 
L/a, where L is the length of the wires. However, Eq. (pi) 
is only an approximation that holds in a limited range of 
T. Such an apparent power-law behavior has been ob- 
served in a number of ID and quasi-lD systems.'^^HISl We 
demonstrated that the RCH provides a more plausible 
explanation of these observations than models based on 
the concept of ID Luttingcr liquid-'^' 

The present paper is devoted to the crossover from the 
low-T RCH to the higher-T Mott law. To our knowl- 



edge this problem has not been studied theoretically al- 
though similar problems have been examined in the con- 
text of VRH transport along the longer dimension of 
three-dimensional (3D) filmJI^and two-dimensional (2D) 
strips .li^ 

We show that depending on space dimension d, the 
Mott law and the RCH represent either two separate, 
competing contributions to the transport or they succeed 
one another via a continuous evolution. In the first case, 
realized in d = 3, the logarithmic derivative d\na/dT 
has a sharp change at the crossover point. Systems of 
dimension d < 3 produce the other type of behavior, 
where a gradual variation of Incr(r) takes place. Our 
analysis is most complete in ID (by which we again mean 
an array of ID systems connected in parallel), where we 
can utilize both numerical and analytical methods. More 
complicated cases of higher dimensions are studied using 
qualitative physical considerations. 

The paper is organized as follows. Sections [H] and |III| 
introduce the model and key relations we use to analyze 
the problem. Sections |IV| and |V] deal with d > 1 and 
d = 1 systems, respectively. Concluding remarks are 
given in Sec. |VI| Details of the derivations are gathered 
in the Appendix. 



II. CONDUCTIVITY FROM CROSSING 
PROBABILITY 



We consider localized states (sites) distributed ran- 
domly with the uniform average density g in energy- 
position space. The x-coordinates of the states belong to 
the interval < x < L. Ifd>l, the system is assumed 
to be infinite along the remaining d—1 coordinates. The 
X = side of the system is the source electrode and the 
X = L side is the drain electrode. We wish to compute the 
conductivity it{T, L) between the source and the drain. 

We adopt the Miller-Abrahams resistor network 
modePSl in which the resistance between any two sites 



is given by 



Rij — -Roe"' 



2r,: 



e., - So 



2T 



, (3) 



where rij is the distance between the sites and Ei is the 
energy of state i measured from the Fermi level. The 
coefficient Rq should have some power-law T-dependence 
determined by the electron-phonon coupling. However, 
we will ignore it and treat Rq as a constant. 

The next approximation is to view our system as an ar- 
ray of hypercubes of length L and cross-section area L'^~^ 
connected in parallel. (For an array of d = 1 wires this 
approximation is exact.) It yields the relation between 
the conductivity and the ensemble-averaged conductance 
G of a single hypercube: 



a{T,L) = L'-''{G) 



(4) 



In order to estimate (G) we use the percolation theory ar- 
gument. We assume that the conductivity of a given hy- 
percube is dominated by an optimal subnetwork, which is 
constructed as follows ."JJ^ Each pair of sites with Uij < u 
is considered a connected bond; otherwise, it is a broken 
bond. Gradually increasing w, one reaches some value Uc 
where for the first time the source and the drain become 
joined by a path of connected bonds. The dominant sub- 
network is obtained increasing m up to u^ + 1 or so. Fur- 
ther increase in u is assumed not to lead to a significant 
growth of the conductivity because the added resistors 
would be shunted by those already present in the circuit. 

For each disorder realization, the critical u is a ran- 
dom number. Its statistical properties are encoded in 
the so-called crossing probability^^HSO] p(^u), which is a 
cumulative distribution function of u. Thus, P{u) gives 
the probability of having a connected path (the span- 
ning cluster) between the source and the drain via the 
bonds Uij < u. (Accordingly, the derivative of P{u) is 
the distribution function of the percolation thresholds in 
a finite-size hypercube, the earliest study of which was 
carried out in Ref. [211) Understanding the geometry of 
the spanning cluster is a key to calculating both P{u) 
and the average conductance (G) . 

In general, the cluster consists of the current-carrying 
backbone and the dead-ends. We can imagine two lim- 
its. The backbone can be made of a single filament, as 
illustrated in Fig. Ilja) and (b), or it can look like a d- 
dimensional network with a correlation length (charac- 
teristic size of the cells) ^ <C i, see Fig. fl|c). In the 
latter case, which is possible only if d > 1, we find 



(G) ^ [L/O 



d~2 



-l„-« 



duP'(u)R^^e 



{L/if-'^R-^ / duP{u)e-'' 



(5) 



The second line in Eq. ([5]) is obtained integrating by 
parts. Substituting it into Eq. (HI), we obtain, with the 
same accuracy. 
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FIG. 1. Evolution of a 2D network with increasing T (see 
Sec. IV I. The network progresses from independent conduct- 



ing strands to an interconnected grid, see main text. 

Similarly, in the case of a single filament of a charac- 
teristic transverse dimension L_l, we can estimate the 
conductivity as follows: 



i?n 



max • 



{LLi7^(u)P(M)e-"} 



(7) 



a - R^^ max {i^-'\u)P{u)e-'') 



(6) 



The brunt of the remaining work is to compute P{u) and 
L_\_ (or ^) as a function of T, L, and d. 



III. TWO MAIN REGIMES 

As discussed in Sec. [l] there are two principal regimes 
of VRH, the Mott law and the RCH. In this section we 
show how they follow from our Eqs. ([5|-([7|. 

A. The Mott law 

To see how Eq. ([s]) leads to the Mott law let us re- 
call that according to the percolation theory, in the limit 
L — >■ oo, the crossing probability P{u) approaches a step- 
function, P{u) — ;> Q{u — Uc)', therefore, 

Incr ~ —Uc ■ (8) 

The integral in Eq. ([5]) is dominated by the interval 

Uc < u < Uc + I ■ (9) 



To calculate the threshold value Uc we proceed as follows. 
From Eq. ([3]) we see that for a given u we need to consider 
only sites with energies ~Tu < e < Tu. These sites have 
the average coordination number 

^(") = pT^^-^ ' (10) 



which scales with u as 



d+l 



C{u) = Pd [u/um) 
UM ^ (2To/T)i/('^+i) 



(11) 
(12) 



where /3i = 1/2, h = 7r/8, and /^a = 7r/20. The per- 
colation threshold is known to be C^ ~ 3 for d — 2,3. 



1/(^+1; 



UM- 



Hence, the threshold value is Uc = {Cc/Pd) 
Substituting this into Eq. (Is]), we recover the Mott law 
[Eq. (fll)] and determine the numerical factor therein to 
be Q = 2CJl3a. 

The geometry of the critical subnetwork in the Mott 
regime is of course more complicatecF^ than what is 
sketched in Fig. fl|c). This subnetwork is a self-similar 
fractal object on the spatial scales tm ^r <^ ^. Here 



rM = auul'^ 
is the typical hopping length. 



i = tmUt^i 



is the correlation length, and 






1/2=4/3 (d = 2) 
1.3 « 7/8 (d = 3) 



(13) 
(14) 



(15a) 
(15b) 



are the percolation theory exponents.'^'^ll Equation ( 14 ) 
for ^ follows from the general formula 



where 



i-^TuVX 



U/Uc 



1 



(16) 



(17) 



is the fractional distance to the percolation point. For 
u = Ur + I [cf. Eq. ([9])], we have e ^ ^Jum, leading to 



Eq. (14) 



The above derivation applies if d > 1. In ID the per- 
colation threshold does not exist in the L — > 00 limit 
because the infinite spanning cluster is invariably bro- 
ken apart by fluctuations. We discuss this special case in 
Sec. El 



B. Rare-chain hopping 



It is easy to see that for ^ given by Eq. (14), the con- 



dition L ^ ^ necessary for validity of the Mott law can 
be satisfied only at high enough temperature. At lower 



T, the network geometry must be different and devia- 
tions from the Mott law should appear. ^'^ For the lowest 
T we expect the RCH regime. Let us now rederive the 
corresponding conduct ivitySQ from our formalism. 

The connectivity of the chain of N sites is determined 
by the probability of forming iV consecutive bonds. Ac- 
cordingly, P(u) ^ \C(u)\^ . The lower bound on N is 
2L/au because the length of each hop is less than or 
equal to 



R = ua/2 . 



(18) 



Since C <C 1, the optimal number N^ of sites must be 
close to this bound; hence. 



InP(ii) 



2L 



In 



(v) 



d+l' 



U <^ Um . 



(19) 



Substituting this into Eq. ([7]), we find the optimal u to 
be, with logarithmic accuracy 



2{d+l)L 



urch 



In (^ 
\urckJ 



1/2 



For the conductivity, we find 
Incr ~ - 




1/2 



(20) 
(21) 

(22) 



in agreement with the previous work.l^ Equation ( |22[ ) rep- 
resents a very slow T-dependence compared to the Mott 
law. If the accessible range of T is limited, as is often 
the case in experiments, such a dependence can be easily 
confused with a power-law, Eq. (pi). In fact, Eq. (22) can 



be written as 

0- - T^* . (23) 

Therefore, the exponent a in Eq. (l2| is essentially the 

opti mal nu mber of sites N^, c>c ln~^"[(T/ro)wRt;g] in the 
chain.lMH 

The remainder of the paper is devoted to analyzing the 
crossover between the Mott law and the RCH. 



IV. 2D AND 3D SYSTEMS 

We make a key observation that at finite L, function 
P{u) is nonvanishing even at u < Uc due to some disor- 
der realizations that percolate "early." This creates an 
exponential tail of P{u), which competes with the fac- 
tor e^" in Eq. (I7|. As a result, the optimal u can be 
pushed below the threshold, m* < Uc, implying that the 
transport is governed by subcritical percolation, e < 0. 

The behavior of the crossing probability near the per- 
colation threshold < — e ^ 1 can be understood quali- 
tatively as follows. Correlation length ^ in Eq. ( 16 ) repre- 



sents the characteristic size of the largest connected clus- 
ters (i.e., the largest clusters among those that are not 



yet exponentially rare^^l). This is different from the geo- 
metrical meaning of ^ in the supercritical regime u > Uc 
where it is the characteristic size of the voids in the net- 
work, see Eq. ^ and Fig.llTc). In the subcritical regime, 
the spanning cluster appears when L/^ independent clus- 
ters of size ^ each join together by chance, forming a 
conducting pathwa y. Th erefore, the crossing probability 
can be estimated agi^Eni P[u) ~ e~^/^. Combined with 



Eqs. (fiel) and (Tlh, this estimate yield^ 



InP(u) ~ - 



L 



This implies 



P'{u) oc {uc — u) 



l^d-l 



(24) 



(25) 



Hence, the behavior of P{u) near Uc depends on whether 
Pd is larger or smaller than unity. In 2D, we have 1^2 > 
1, and so the derivative vanishes. In 3D, the opposite 
inequality I'a < 1 holds, and so P'{uc) diverges. This 
dichotomy is the reason for the different manner in which 
the Mott law transitions to the RCH regime in the two 
cases. 

Consider the 2D case first. Per Eq. (|6| the conductiv- 
ity is determined by the maximum of e^"P{u). Using 
Eq. (24), we find that it is reached at u^. such that 



u,-u, = {Tp/Tf 



T, 



where 



P- 



Vd + l 



RCH 

7 
3' 



<T <Tf 



(26) 



{d+l){vd-l) 
and the characteristic temperatures Trch and Tp are 



(27) 



T, 



Tp^ro(-) 

-a\('i+i)/2 



/a\9/7 



1) /a\9/7 



RCH 



-'& 



To 



L 



3/2 



(28) 
(29) 



Tp is temperature below which the T-dependence of the 
conductivity and the network geometry start to deviate 
from what is found in an infinite sample, cf. Sec. |III A| 
and Fig.flVc). The Mott law acuires the correction term 
as follows: 



In cr 






(Uc 



Formula (14) for ^ is replaced by ^ ^ L/{uc 



(30) 

t*) < 

L and the VRH network now consists of well separated 
filaments, see Fig. [l|b). 

As T decreases further and reaches Trch, the optimal 
u becomes equal to urch [Eq. ( 21 )]. At this point the dis- 
tance from the critical point u^ — u^ becomes comparable 
to Uc itself. The cluster size ^ shrinks down to the ele- 
mentary hopping length i?, which signifies the crossover 
to RCH, see Sec. |IITb| 

In 3D, the optimal u* has a different T-dependence. 
Due to positive concavity (the second derivative) of the 
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FIG. 2. Crossing probability P{u). The qualitative difference 
of the (a) 2D and (b) 3D cases in the vicinity of Uc is apparent. 
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FIG. 3. Conductivity as a function of temperature. The tilted 
dashed line shows the extrapolation of the Mott law past its 
region of validity for contrast, (a) In 2D, there is an inter- 
mediate regime smoothly interpolating between the Mott law 
and the RCH. (b) In 3D, the transition between the Mott and 
RCH regimes is abrupt. 



InP(M) curve near Uc, see Fig. ^ the maximum of 
e~"P{u) stays at Uc at all T > Trch- Temperature 
scale Tp plays no role and Eq. ([26]) does not apply. The 
Mott law continues to be valid. However, as soon as T 
drops below Trch, parameter u* gets suddenly reduced 
by some numerical factor and at still lower T we have 
Eq. (20). In other words, in 3D the transition from the 
Mott law to RCH is abrupt, see Fig. [3t 



V. QUASI-ID SYSTEMS 

A. Toy model 

Since the percolation threshold does not exist in ID, 
function P{u) for high u has a somewhat different form 
compared to d > 1. The peculiarities of the ID case are 



well illustrated by the toy model in which all the site 
energies are equal to zero. The connected sites are those 
that are separated by distance less than R [Eq. (18)]. 
The crossing probability can be computed exactly: 



^ = 1+ E ^-^^^''w-k)] 



fe-1 



(k-kp + lp), (31) 



l<k<l 



where I — L/R, p — gR, and the density of states g is 
redefined to be simply the density in coordinate space. 
For large I this expression can be approximated by 



lnP^iso-p)l, 
where sq is the real root of the equation 



m 



l-e- 



1 
P 



cf. Appendix [A} For large p, this yields 

\nP~-lpe-P , p>l. 
For small p, we have sq ^ Inp, which gives 

InP ~ -Hn- , p < 1 . 
P 



(32) 



(33) 



(34) 



(35) 



These results have a simple interpretation. For high p, 
the sites are very dense. However, empty segments of 
length larger or equal to R may also appear by chance. 
Deviations of P from unity are due to these rare disrup- 
tions — "breaks." The probability of having no sites in 
a segment of length R is 



Ptoy 



The average distance l/(ptoy5) between such breaks de- 
fines the average length ^ of connected clusters. Using 
InP = -L/^, as in Eq. ptl), we recover Eq. ([34]). On the 
other hand, for small p, the sites are very dilute, so we 
have an analog of the RCH regime. We expect P ^ C^, 
similar to Sec. |IIIB| Approximating the number of hops 
N by I and the average coordination number C by p, we 
arrive at Eq. ( 35 ). (We assume that N is large and ignore 



its fractional part in this heuristic derivation.) Note that 
Eqs. (34) and (35) match at Inp w —1, where the cluster 
size shrinks to the elementary size R. In this respect, our 
toy model is similar to the 2D case (Sec. IV). 

Next, we show that the having random energies in 
addition to random x-coordinates does not qualitatively 
change this physical picture. 



B. Analytical results for ID 
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FIG. 4. VRH network in ID. (a) At low T, the network is 
made of a few approximately equidistant hops. The resistance 
of each link e" is much smaller than the typical one e"*^ . This 
structure is common to all d, cf. Fig. [If a), (b) At high T, most 
of the resistors in the network (dashed lines) have typical 
values and form clusters. However, unlike the d > 1 case 
[Fig. [Ub)], the clusters are connected together by u ^ um 
"breaks" that straddle regions empty of the hopping sites. 
The solid line shows the path of the least total resistance. 



the limits of either high or low T where a is dominated 
by rare events. Below we rederive these results in a uni- 
fied manner, which enables us to elucidate the crossover 
between them. 

Recall that each link of the hopping network is char- 
acterized by a dimensionless number Uij. As in the toy 
model, we can talk about "breaks," by which we mean 
links with Uij much larger than the typical value um- 
At hi gh T the transport is dominated by these rare 
breaks ,'^^'2lH2ll see Fig.[4]jb). The breaks straddle regions 
of size R = ua/2 and 2uT in x and e directions, respec- 
tively, empty of hopping sites. The probability of having 
(36) a break with Uij > u can be computed considering how 



the area and the perimeter of such regions scale with u. 
The result iJ^HSl 



p{u) 



exp 



u 

"if 



2B- 



um 



(37) 



whereas B w 0.9. The average size of the connected clus- 
ter is therefore ^ = vm /p{u) and the crossing probability 
is, similar to Eqs. ([24|) and (|34|), 



InP(u) 



2L 



exp 



2- +25 



(38) 



For such P{u), the maximum of the right-hand side of 
Eq. ([?! is reached at 



u* ~ ujv/ln^ 



-'RCH 



"M 



Um < wrch 



(39) 



The model of ID VRH in which both the coordinates 
X and energies e of the sites are random has been studied 
extensively in prior literature. It has been shown that a 
number of rigorous analytical results can be obtained in 



Accordingly, the conductivity a is given bjSE 

1/2 



lna(r)~- 



270^ /LT 
T la To 



(40) 




tion ( 44 1 is also in agreement with the general RCH for- 
mula (19), and so it leads to Eq. (22 1 within the leading 
logarithmic approximation. In principle, Eq. (42 1 goes 



to the next order in this approximation. In practice, the 
logarithms involved are never truly large. This moti- 
vates us to additionally examine the problem by numer- 
ical simulations, which are described in the last part of 
this Section. 



FIG. 5. Transport regimes in the quasi-lD case. The lowest 
T region (to the r ight of the second dot on the curve) is the 
RCH regime, Eq. (221. The intermediate T-regime (between 
the two dots) is described by Eq. ( |40[ ). It transitions smoothly 
into activated transport at high T. Near the inflection point 
(the leftmost dot) the Mott law, Eq. ([l|, is realized. 



It was shown in Ref. [57|that this equation has the upper 
limit of validity, T = ro/(4ci), where ci is defined by 



Cl 



21nu2 



RCH 



2 In 



2L 
a 



(41) 



At still higher T the conductivity is described by the ac- 
tivation law,25 27 incr(r) ~ -To/(2T). This imphes that 
Mott law is, strictly speaking, invalid in ID. However, if 
the accessible range of temperature is narrow, the devi- 
ations from the Mott law may not be readily apparent. 
These deviations are the smallest precisely near the point 
of crossover to the activation law because the dependence 
of Incr(r) on um oc T~^''^ must have an infiection point 
there. In this sense, we can argue that the Mott law 
is confined to a narrow vicinity of T = To/(4ci), while 
Eq. ( 40 1 describes an intermediate subcritical percolation 
regime, similar to Eq. (30) in 2D. We may expect that 



this intermediate regime should cross over to the RCH 
when the cluster size ^ becomes of the order of the aver- 
age hopping length tm ■ Indeed, this criterion is satisfied 
at Um ^ urch, see Fig. [5] 

For u < urch, the behavior of P{u) is determined by 
the RCH, see Fig. [Ij^a). Adopting the method of Ref. E] 
to the ID case, we obtain the result conceptually similar 
to Eq. (laSl): 



InP(u) 



2L 



co(u) 



(42) 



Function cq{u) in this equation is defined in Appendix [b| 
For u ^ Um, it has the asymptotic form 



co(u) 



1 



In 



'M 



SO that the crossing probability is given by 



InP(M) 



2L 



In 



'■M 



(43) 



(44) 



As in our toy model, P is the product of the number 
of hops N = 2L/{ua) and a logarithmic factor. Equa- 



C. Numerical results for ID 

To compute P{u) numerically we proceed as follows. 
First, we generate an ensemble of ID systems for each 
L and T of interest. For every realization we find Uc 
— the minimum value of u at which the path traversing 
the system is formed. This Uc is a random number be- 
tween and 2L/a. To find it we start from some initial 
guess ~ L/a and then fine-tune it by iterative bisection. 
The connectivity of the network for a given u is deter- 
mined by t he a lgorithm similar to one used in our pre- 
vious work.lSESl Finally, the histogram of u^ gives us the 
crossing probability P{u). 

In Fig. p^ we show that when our results for P{u) are 
plotted as\r M / L) InP(u) versus u/um they collapse on 
a common master curve. The deviation from the master 
curve at high u come from the fact that the number of 
hops A^ ^ 2 is no longer large. The deviations at low u 
are due to poor statistics of these low probability events. 

We considered two analytical approximations of our 
numerical results. First, aiming to utilize the analogy 
to the exactly solvable toy model, we fitted the master 
curve in terms of function P defined by Eq. (31 ). The fit 
(shown by the longer curve in Fig. pi) was generated by 




FIG. 6. (Color online) Numerical results for P{u) in ID. 
Difi'erent colors represent T ranging from 5 x 10~^ to 5 x 10~^ 
for L = 100a (circles) and from 10"" to 10"^ for L = 30a 
(squares), with the temperature unit being such that To = 
3/4. The solid lines are analytic approximations, see main 
text. 
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FIG. 7. (a) ID conductivity vs. um obtained numerically 
(see main text). The tilted dashed line and the two dots have 
the same meaning as in Fig. [5] The parabola on the left 
of the dashed line represents the activation regime, (b) The 
same data plotted on a double logarithmic scale in order to 
demonstrate that it can resemble a power-law (pi at low T 
(to the right of the dot). 



setting 



1 



1 + 2B2 



(45) 



'■M 



designed to produce p ~ ptoy at u ^ um, cf. Eqs. (36) 
and (|37|). 

The second approximation is the parameter-free 
Eq. (42), which simultaneously represents the leading 



asymptotic result for P{u) at low u and the strict lower 
bound at any m, see Appendix [B] It is shown by the 
shorter solid curve in Fig. [6] While our numerical re- 
sults converge to this curve at lowest m, we were unable 
to achieve an accurate match because of prohibitive com- 
putational cost needed to acquire statistics in this region. 
Using our first fitting formula for P{u) and Eq. ([5]), we 
also computed the conductivity at sufficiently large um > 
12 where this approach should be valid. The results are 
shown in Fig.JTJtogether with the activated dependence^ZI 
at Um < 6 and the dashed line through the surmised 
inflection point at um = y/8ci ~ 8. The obtained graph 
suggests a smooth transition from the exponential to the 
stretched-exponential to the power-law-like behavi or o f 
(t(T), consistent with Fig. ^ and our previous work.!^ 



demonstrated that both stretched exponential and quasi- 
power-law dependences can arise in such systems. The 
significance of our results in elucidating the transition 
between these different laws. 

Prior theoretical work that dealt with hopping conduc- 
tivity across a thin film took the n umber of hops N as the 
principal variable of the problem.'^S^Zl This allows one 
to compute the conductivity in the RCH regime beyond 
the leading logarithmic approximation given in Sec. |IIIlJ| 
However, this approach precludes one from making con- 
nection with the percolation theory and the Mott law. 
Our method, which employs concepts of subcritical per- 
colation theory and takes u as the main variable, allows 
us to do so. In the RCH regime, the computation of the 
optimal u = u^, implicitly optimizes the number of hops 
in a path. In the Mott regime, it plays the role of the 
percolation parameter directly. 

Our results provide a natural explanation for the be- 
havior observed in quasi-lD systems with a large number 
of channels .liSHill However, individual wireP^and similar 
systems, like graphene nanoribbons, exhibit significant 
mesoscopic fluctuations. Further effort is required to in- 
clude these fluctuations into account. 

Examples of 2D systems that can be studied using 
our method include GaAs devices^ and bilayer graphene 
p-n junctions.^ Finally, in 3D our approach describes 
thin films. ^ "^ ^^ Reference ,3; indeed reported the crossover 
from Mott to a power-law-like dependence with decreas- 
ing temperature. 

In an interesting recent experimenPS the transport 
across disordered films of Cu-phthalocyanine has been 
measured over a broad range of temperature T and 
thickness L. A sharp change from strong to weak T- 
dependence below some temperature in the thinnest films 
studied therein agrees qualitatively with our theory. The 
dependence on L is more complicated to analyze because 
of possible systematic variation of the film morphology 
and doping level with thickness. 

Extensions of our work may include the study of 
Coulomb and spin-related effects, I^SHSII energy-dependent 
(e.g., ex ponen tial) density of states, and non-Ohmic 
transport .l^Eni 

This work is supported by the University of California 
Office of President (UCOP) Program on Carbon Nanos- 
tructures and ACS of UCSD. 



Appendix A: Solution of the ID toy model 

Within our toy model, the (unnormalized) probability 
of connecting the source and drain via A^— 1 intermediate 
sites is given by 



VI. DISCUSSION AND CONCLUSION 

In this paper we have investigated the dependence of 
VRH conductivity on temperature in systems of finite 
thickness or length in all physical dimensions. We have 



^ N N 

where Xi is the length of the ith hop and g is the density. 
The normalized probability is obtained by summing over 



N and multiplying by the factor e~^^ known from the 
Poisson distribution. In order to compute each term we 
first take its Laplace transform over L. This decouples 
the integrals and yields the closed-form expression 



Zn_i{R,s) =5 



N-l 



1-e 



-sR 



N 



(A2) 



Next, we perform the inverse Laplace transform and sum 
over N. The result is 



N 



^w = EE "''":!°if'""""' e(^-^") 



N=l fc=0 



kl{N-k)l 



(A3) 

In terms of the notations I — L/R and p — gR introduced 
in Sec. IV) the crossing probability becomes 



P = 



oo N 

EE 

^=1 fc=0 



N-l„-i: 



N\ i -l)^[p{l-k)r -^e 
{N-l)\ 
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which can be shown to lead to Eq. ( 31 ) by straightforward 
algebra. 



One can also proceed by summing Eq. ( A2 ) for all N 
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The inverse Laplace transform can be written, after the 
change of variable s —> s/R, as 



-pi 
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2m 
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(A6) 



The integration contour C is shown in Fig. [8] 
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FIG. 8. The contour used in Eq. (|A6|. For p < 1, all the 
poles Sj are in the left half-plane. When p becomes sufficiently 
large, the pole so, which has the largest real part becomes 
positive and approaches p. 



Applying the residue theorem, we get 



p- E 



p{l- p + s^) 



^{Sj~p)l 



(A7) 



where Sj are the roots of Eq. ( 33 1 . For large I the real 
root So dominates; therefore. 
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So 



p(l-p + so) 



^(so-p)l 



(A8) 



Equation (32 1 quoted in Sec. Iv] was obtained from here 
by dropping the pre-exponential factor. 



Appendix B: Lower bound on the crossing 
probability 



For any d we define the function 
•^ \ fc=i / i=i 

N N-l 

dV ^^dxd'^-^ri Jl gdsj 



u-Ui,j_i), (Bl) 
(B2) 



i=l 



Unlike Zf^^i in the toy model [Eq. (All] here we take 
into account random site energies —Tu < Ej < Tu. The 
energies of the electrode sites, eg and ejq, are both equal 
to zero to minimize the resistance. The transverse dis- 
placements r~" that exist in d > 1 case are unrestricted 
except the first and the last one, r^ = r;^ = 0. The hop 
lengths Xi in the x-directions are assumed to be positive. 

The sum of Z^-i over all N yields the unnormalized 
probability of finding a path with Ui^i-i < u for all links 
i provided only the hops between sites that are nearest- 
neighbors in X are allowed. Hence, this quantity gives 
the strict lower bound on the crossing probability P{u). 
The bound becomes sharp at u <C um where the average 
number of nearest neighbors is parametrically small. 

Similar to Appendix|Aj we proceed to take the Laplace 
transform of Zf^^i with respect to L: 
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Zn-i= d^W^ 



-L\2 I ^2 ^~sx 
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(B3) 
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Henceforth we focus on the ID case. Making two changes 
of variables, C = e/{uT) and c — as/2, and integrating 
over Xi, we obtain 



N 



Zn-i^{^) (guT)''-' fr,MO) , (B5) 



where /o(C) = ^(C)j function k{() is defined by Eq. (33 1, 
and functions fj (C) with j > obey the recursive integral 
equations 



n-ce(-c) 
/,(C)= / df^AiC,rj)f,^,{7j), (B6) 



i+ce(C) 

A 



A(C,7?) = l-e 



(B7) 



A(C.)-ll- '^' + '^':'^'"' h^- (B8) 



Similar equations appeared previously in Ref. ini a near- 
literal copy of which is available online as Ref. |40l Fol- 
lowing this work, for A^ ^ 1 we expect fj^_i{Q ~ 
aK'^~^tp{C,), where k = >c{uc) is the largest eigenvalue 



of the integral operator in Eq. ( B6 ) and a is the overlap 



of /o(C) with the corresponding eigenfunction ipiO- This 
approximation enables us to derive the leading asymp- 
totic behavior of the inverse Laplace transform of the 
geometric series: 



-iC30+7 



ZTTl I a ^-^ 
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— ioo+7 



Representing this integral by a sum over the poles of the 
integrand, we see that the dominant pole cq is the real- 



valued solution of the equation 

MCo 



^J\/ 



(BIO) 



while the crossing probability is given by Eq. (42). For 



u <^ um the product ucq proves to be large and negative. 
In this limit we can reduce our integral equation to that 
studied in Ref. ^, and obtain the eigenvalue equation in 
the form 



o\uco\ 



\UCq\ 



-kt , >ct ~ 1.18 . 



(Bll) 



Solving this equation in the leading logarithmic approx- 
imation yields Eq. (43). In order to produce the second 
analytic fit (the shorter solid line) in Fig. [6] we solved 
both the integral eigenvalue problem and Eq. ( |B10 ) nu- 
merically. Hence, we are quite confident that what is 
shown in that Figure is the strict lower bound on P{u). 
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